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Abstract 

We solve exactly a class of Frenkel-Kontorova models with piecewise parabolic 
potential, which has d sub-wells in a period. With careful analysis, we show 
that the phase diagram of the minimum enthalpy configurations exhibits the 
structure of a complete d-dimensional devil's staircase. The winding num- 
ber of a minimum enthalpy configuration is locked to rational values, while 
the fraction of atoms in each sub-well is locked to values which are sub- 
commensurable with the winding number. 
PACS number(s):03.20.+i, 05.45.+b, 64.60.Ak 
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Periodic modulated structures are quite common in both condensed-matter and dy- 
namical systems. As suitable parameters are varied, such structures may go through 
commensurate-incommensurate phase transitions. In particular, when a system is in the 
highly nonlinear regime, there is generally a tendency for the periodicity to lock into values 
which are commensurable with the lattice constant The Frenkel-Kontorova (FK) model 
is just such an example, which describes a one- dimensional chain of coupled atoms subject 
to a periodic potential V: 

H{{u n }) = \ ("n+1 - u ri f + XV (u n ) - a(u n+l - u n ) . (1) 

Here, u n is the position of the nth atom, A is a measure of the nonlinearity. For simplicity, 
suitable length scale is chosen so that the periodicity of V{u) is one. 

From the variational principle, we see that a stationary configuration {u n } must satisfy 
the following equation: 

u n+ i - 2u n + u n -i = XV'(u n ). (2) 

It is well-known that the above equation can be rewritten as coupled first order difference 
equations, which provides a link between the stationary configurations in an FK model and 
the trajectories in a two-dimensional area-preserving map ||. 

According to Aubry [§,|J, a minimum-energy configuration is one in which H cannot 
be decreased by alteration of a finite number of u n 's. For these configurations, there is a 
well-defined winding number 



u N -u_ N < 

00 = hm -— — — , (3) 

N,N'^oo N + N' 



which is the mean distance between neighboring atoms. A ground state is a recurrent 
minimum-energy configuration, which can be depicted by HQ 

u n = f u {nu; + a). (4) 

Here, a is an arbitrary phase variable; the hull function f u (x) is increasing with x and 
satisfies 



Ux + l) = Ux) + l. (5) 
When fu>(%) is discontinuous, the positions of the atoms is described by either fu(x + ) or 

fu(x~). 

There are two distinct types of phase transitions in the FK models with smooth potential. 
First by varying a with A fixed, one finds a commensurate-incommensurate phase transition. 
For a given a, there is a uj(ct) which specifies the minimum enthalpy state and it takes on 
rational value on an infinite set of plateaus. Such structures are called incomplete or complete 
devil's staircases depending on whether or not there are incommensurate phases with non- 
zero measure between the commensurate ones. There exists a critical value of A, 
above which the devil's staircase becomes complete. On the other hand, if u is fixed while 
A is varied, one will see that for the incommensurate ground states there is a non-vanishing 
Peierls-Nabarro barrier above a critical value of A It is called the transition by breaking 
of analyticity (TBA), since the most evident indication of this transition is the change in the 
analyticity of the hull function [[|. In terms of the two-dimensional area-preserving map, 
what happens to an incommensurate ground state is that it changes from a configuration 
represented by a KAM torus to one which is represented by a Cantor- Aubry-Mather (CAM) 
set §. 

The first exactly solved FK model was studied by Aubry using a piecewise parabolic 
potential [@,|5]|§. The potential is the simplest non-linear one in that the non-linearity only 
occurs at one point in a period. It turns out that in this model both the critical values of 
A are zero for the transitions described above. In other words, the devil's staircase in this 
system is always complete and the incommensurate configuration is always pinned unless 
A = 0. Simple as the model may seem, it provides us a clear picture of systems in the pinned 
phase. 

In this paper, we will solve exactly an extension of the Aubry model, when the potential 
has d sub- wells in a period. Such substructure may occur naturally in surface reconstruction 
or be realized in a superlattice. This model was first proposed by Griffiths, et al [0. Several 
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interesting new phenomena such as solitons in the incommensurate phase, discontinuous 
Cantorus-Cantorus phase transitions, etc were found in the d = 2 case. As we shall see, the 
phase diagrams in the higher dimensional case are much more complicated 
When there are d pieces of parabolas in a period, the potential is given by 

V{u) = min {Vi(u)}, (6) 

l<t<a 

with 

Vi(u) = ~ {u - hi - Int[w]) 2 + hi. (7) 

Here, Int[w] is the integer part of it; hi and fej's are generally independent parameters. 
Requiring V(u) to be continuous, we find that the positions of the kinks are given by 

h + bi+i h i+ i — hi 

= —r- + ^r^7- < 8) 

with hd+i = 1 + hi and b^+i — 1 + bi. p^s must satisfy 

< Pi < . . . < p d < I. (9) 

For convenience, we will choose pd = 1, which leads to the constraint hi + b\j1 = hj + (1 — 
bd) 2 /2. An example for d = 3 is shown in Fig. 1. 
The force-balance equation is given by 

u n+ i + w n _i - (2 + \)u n = -A (Int[w n ] + , Pi- X < u n - Int[u„] < Pj. (10) 

Introducing the hull function f w (x) [|TT| and making use of the fact that it is increasing with 
x M, we can rewrite the above equation as 



f u (x + + f u (x -u) - (2 + X)f u (x) = -X lh + J2(b i+1 - bi)lnt[x + 1 - A]| , (11) 

where fa = Y^j=\ v y v % designates the fraction of atoms that lies in the i-th sub- well and 
Z)f=i v i — 1- F° r consistency, /3j must satisfy the following condition: 

upn<^<upt)- (12) 



It is easy to verify that the solution to Eq. (|TTD is given by 

oo ( d ~\ 

Ux) = d E e " Hx &i + E(^+i-^)Int[n^ + x + l-A] • (13) 



Here e" "* = {l + A/2 - y / A(l + A/4) }, and d = tanh(x/2). It is obvious from the above 
result that around the i-th kink at Pi there will be a region of length (b i+ i — ) d , in which 
no atom is allowed. 

Once we know the hull function, it is straightforward to find the average enthalpy per 
atom: 

(u 2 - 2au) A 



ij(u, Ui,..., v d ) 
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~ E Uk-b^l + k-b^vX + xj^hi 

Z Ki<j<d L J i=l 



A d 



d 



I E e "' n|x E (ku - Wi+i " &i) ^(na; + & 4 ) . (14) 

n=— oo lij=l J 

Here, S'(x) = (x — Int[x]) 2 — (x — Int[x]) , and f3j t i = /3j — The hull function and the 
average energy per atom have been obtained independently by Griffiths, et al in a different 
form |TJ. 

Since rj(uj, V\, . . . , vj) is a convex non-differentiable function, the minimum enthalpy con- 
figuration for a given set of (a, hi, . . . , ha, £>i, . . . , &d) must be determined by the condition 
that the directional derivatives are non-negative in all directions. After some work, the 
condition simplifies to 

A d 



cos 9 < —a + 



oo , d—1 s I 

E ne-l^(E A& i A6,Int[no; + ^ J + ^]) 

n=— oo J,fc=0 J 



+ A £ cosflJ A/i* + 4 E e" |n|x (E E Aftj-Aifclnttnw + (5 k)j + a»J ) ) > 0. (15) 

* — 1 ' n=-oo 0=0 k=i 



k j 

Here, Ahi = (/i, — /^), A5,- = (&,+i — &,•), = e(ncos6 ) + £ cos6*; — £ cosfy). cost's 

z=i z=i 

satisfy J2 cos 2 0j = 1, and < &i < it. 5% j is just a bookkeeping device to remind one 

i=0 

of which limit to take in the function Int[x] at integer values of x. It is obvious from the 
above result that the phase diagram only depends on (Ah\, . . . , Ahj-i, A&i, . . . , Afa^-i), as 
it should be. 



To illustrate how the phase diagrams look like, let's first concentrate on the case of d = 2. 
For a given set of (u, Ahi, Abi), there is only one v\ that gives rise to the minimum energy 
state, and it can be determined by setting 8q = ±7r/2 in the minimum enthalpy condition, 
Eq. (f0|). The resulting v\ is locked to values sub-commensurate to to, i.e. there exists some 
integer n such that nuo + V\ G Z (the set of all integers). In Fig. 2, we show the phase 
diagram for u = (\^5 — l)/2. Each vertical cross section will give a devil's staircase in the 



Ahi-Ux plane. A special case for A&i = 1/2 has been shown in Fig. 2 of [10| . 

The phase diagram of the minimum enthalpy states in this case is three dimensional. For 
simplicity, let's look at the constant Abi and Ahi slices, respectively. In the former case, 
the domains of stability fall into three categories (see Fig. 3): 

i. ) n u, {riiuj + Ui) G" Z, for any n , rii G Z (uj irrational): 

The domains are discrete points. 

ii. ) There exists n G Z such that either nu or {nu + U\) G Z. 

1. u = p/q ||13|| . (nuo + v\) G^ Z for any n G Z: 

The domains are horizontal lines with length (2A6? - 2A 6l + IjAdpge^ ^ 

(1 — e qx ) 

2. irrational, (ncj + Vi) G Z for some integer n : 

The domains are lines with slope — l/(An) and length 
A6i(l - A6 x )Vl + A 2 n 2 (d e~ |n|x ) . 

iii. ) ui = p/q, (su + Ui) G Z for some integer s G [0, g — 1] : 

The domains are infinite-sided polygons with sides given by the lines described in (ii) 
and vertices by the points in (i). 

Summing the area of all these polygons, we verify that they form a two-dimensional 
analog of a complete devil's staircase in the a — Ah\ plane. The phase diagram we obtained 



here is very similar to that obtained by Aubry, et al in a different model |TJ]] and that 
obtained numerically by P. Delaly [15[ . 



For constant Ah\, the phase diagram is very similar to the previous one, except for the 
following two aspects (see Fig. 4): 

a. ) One can easily see from Fig. 2 that only some of the vis satisfying the sub- 

commensurability condition are allowed; 

b. ) the lines in (ii.2) are parabolas instead. 

It is straightforward to carry out a similar analysis for the phase diagrams in higher 
dimensions, i.e. d > 3. Let us introduce a = a + 8a and Ahj = Ahi + 8hi. Here, 

\J oo , d— 1 >. 

a = J2 ne-^\ ]T A6 3 -A6 fc Int[nc; + ^ fcJ ] , 

n=-oo j,k=0 * 

oo , t— 1 d— 1 n 

Ahi = d e' Hx EE^A^ + ^il > ( 16 ) 

n=— oo jf=0 fc=i 

with Int[x] = (Int[x + ] + Int[x~])/2. Eq. (|T5|) becomes 

"■ V <^ E e " Hx | E ^Ab k Z{nw + p k j)\Vl r n\ . (17) 

Here, ^(x) takes the value 1 or depending on whether x is an integer or not. 
u = (cos^q, cos#i, . . . , cos6*d_i), an arbitrary unit vector, and 

v = (8a, —X8hi, . . . , — \8hd-i). Vjfj = (n, 0, . . . , 0, 1, . . . , 1, 0, . . . , 0), where the value 1 is 
assigned to the [min(j, k) + 2]-th to [max(j, k) + l]-th components. It can be seen that for 
a domain of stability to have non-zero measure, the corresponding uj must be rational and 
all the Ui's must be sub-commensurate to to. Furthermore, these domains are enclosed by 
infinitely many hyper-planes, where either u is irrational while all the sub-commensurability 
conditions are satisfied or u is rational while exactly one of the sub-commensurability condi- 
tions is violated. The intersections of these hyper-planes in turn give the boundaries of even 
lower dimensions, where two or more of the commensurability and sub-commensurability 
conditions fail. They correspond to multi-critical points in the phase diagram. 

In this letter, we solve exactly a class of extended solvable FK models, whose potential 
has d sub-well in a period. Both the hull function and the enthalpy per atom for an arbitrary 
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d are given explicitly. Since the enthalpy per atom is non-analytic, the global minimum is 
found by requiring all the directional derivatives to be non-negative. Using this condition, 
we are able to determine analytically the complete phase diagram. For given Ab^s, the phase 
diagram is an extension of the devil's staircase to d dimensions as long as A > 0. It will be 
interesting to see if the conclusion is still valid in the d — > oo limit. Detailed discussion of 
the phase diagram and the associated physics will be given in a subsequent paper. 
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Fig. 2 Phase diagram in the Ab±-Ahi plane for A = 0.1, and uo = (y/E — l)/2. The values 
of n are inset, which is related to v x through v\ — 1 + Int[no;] — nuj. The two heavy lines 
give the upper and lower bound of Ahi. 



12 




Fig. 3 Domains of stability (infinite-sided polygons) for (A, A&i) = (0.5, 0.25) are shown 
up to q = 10. The values of (u>, V\) are inset. 
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Fig. 4 Domains of stability for (A, Ahi) = (0.1,4.5 x 10 2 ) are shown up to q = 13. The 
values of cu and v are put on the top and inset, respectively. 
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